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Electronic and optical properties of potential solar absorber Cii3PSe4 
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We theoretically investigate the electronic and optical properties of semiconductor Cu 3 PSe4. We also report 
diffuse reflectance spectroscopy measurements for Cu3PSe4 and CU3PS4, which indicate a band gap of 1.40 
eV for the former. Hybrid functional calculations agree well with this value, and reveal that the band gap 
is direct. Calculations yield an optical absorption spectrum similar to GaAs in the visible region, with 
a > 5 x 10 4 cm -1 for A < 630 nm. We conclude that the optical properties of Cu3PSe4 are within the desired 
range for a photovoltaic solar absorber material. 

PACS numbers: 71.20.Nr, 78.40.Fy, 78.20.Bh 

The need for low cost solar photovoltaic cells continues 
to drive a global effort to find, evaluate, and refine ma- 
terials that offer both cost effective and highly efficient 
solutions. Within the search among the inorganic semi- 
conductors, one of the most exciting families has been the 
ternary and quaternary copper chalcogcnidcs consisting 
of CuCh4 and ACh4 tetrahedra, where Ch is one of (S, Sc, 
Te) and A represents one or more other elements. Exam- 
ples include the commercially widespread Cu(In,Ga)Se 2 
(CIGS), as well as Cu2ZnSnS4, which has the benefit 
of consisting of low-toxicity, readily available elements. 
One material in this family which has not received much 
attention in this regard is CU3PSC4. The structure of 
Cu3PSe4 has been determined 1,2 , and a photoelectro- 
chemical analysis 3 of isostructural compounds CU3PS4 
and CusPSsSe has indicated indirect bandgaps of 2.38 
cV and 2.06 eV respectively. However the bandgap of 
Cii3PSe4 itself has not been reported previously. In 
this Letter we examine this material, theoretically in the 
scope of low temperature and zero defects, and experi- 
mentally in the scope of room temperature measurements 
on powder samples. 

In our synthesis of Cu3PSe4 and CU3PS4, elemental 
powders of Cu (Cerac, 99.5%), P (Alpha Aesar 99%), 
Se (Cerac, 99.6%), and S (Cerac, 99.999%), are ground 
and sealed in evacuated fused silica tubes. Polycrys- 
talline powders are obtained via solid state reaction 
at 480 °C (Cu 3 PSe 4 ) and 600 °C (Cu 3 PS 4 ) for 24 h. 
The XRD patterns, collected on a Rigaku Ultima IV, 
of the reacted samples are similar to the ISCD card 
#095412 for Cu 3 PSc 4 and #412240 for Cu 3 PS 4 . Hole 
majority carrier type is confirmed by positive thermo- 
electric voltage measurements. The optical bandgaps 
are determined from diffuse reflectance measurements 
(see Figure 1) collected with a spectrometer equipped 
with an Ocean Optics HR4000 UV-VIS detector and a 
balanced deuterium/tungsten halogen source (DH-2000- 
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BAL). The Kubelka-Munk model is employed to de- 
termine the optical absorption edges of 1.40 cV and 
2.38 eV for Cu 3 PSe 4 and Cu 3 PS 4 , respectively. The 
Cu3PSe4 bandgap is within the desired range for solar 
photovoltaic devices (1.0 eV to 1.6 eV), while the Cu3PS 4 
bandgap agrees precisely with Ref. [3]. 

Calculations based on density functional theory (DFT) 
are performed using the projector augmented wave 
(PAW) method 4 as implemented in the plane wave code 
VASP 5 . For the exchange-correlation functional we use 
the generalized gradient approximation (GGA) in the 
PW91 paramctrization 6 for accurate total energy calcula- 
tions. To avoid the bandgap underestimation common to 
standard DFT calculations 7 , we use the Heyd-Scuseria- 
Ernzerhof 8 (HSE) hybrid functional, which combines 
Hartree-Fock (HF) exchange with GGA exchange and 
includes an empirical shielding of the HF exchange. The 
HSE functional has yielded respectable estimates for 
semiconductor bandgaps 9 , although it is generally not as 
reliable as methods based on the GW approximation. 10 
We will see below that HSE performs particularly well for 
calculations on Cu3PSe4. For brevity and the expressed 
application interest, we present theoretical analysis for 
Cu3PSe4 but not Cu3PS 4 . Both materials crystallize in 
the wurtzite-based enargite structure with a simple or- 
thorhombic unit cell with 16 atoms admitting the space 
group Pmn2i (FIG. 1). For electronic and optical cal- 
culations, we use experimentally determined lattice and 
atomic parameters. 2 Ionic and lattice relaxations are dis- 
cussed below. In all calculations, we carefully consider 
the convergence of the result with respect to plane wave 
basis cut-off energy, fc-point density, and time saving ap- 
proximations. 

The HSE determined band structure is shown in 
FIG. 2. The HSE calculations used a T-centered 8x8x8 
fc-point grid with 2x down-sampling for one of the 
fc-point summations in the evaluation of the Hartree- 
Fock exchange potential. 11 The band structure plot it- 
self (FIG. 2) is formed by interpolation of the 8 3 fc-point 
grid using a Fourier method. 12 The band structure most 
importantly shows a direct bandgap of 1.38 eV at the T 
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point, in excellent agreement with the experimental value 
of 1.40 eV. The band structure also shows low lying val- 
leys (<~ 0.1 eV) in both the valence band (at Y) and the 
conduction band (along T-Z). We note that an HSE cal- 
culation for CU3PS4 at experimental parameters 13 gives 
these two valleys as being the valence band maximum 
(VBM) and conduction band minimum (CBM), with the- 
oretical and experimental 3 indirect bandgaps of 2.55 eV 
and 2.38 eV respectively. We calculate electron and hole 
effective masses of 0.17 and 0.60 times the electron mass, 
respectively. Spin-orbit coupling was verified to be small 
in a single GGA calculation which showed VBM split- 
ting and a bandgap decrease of 0.030 eV relative to the 
spinless GGA bandgap of 0.29 eV. 

An analysis of the partial density of states indicates 
that the valence bands from the VBM to about —5.5 eV 
below (see FIG. 2) have similar orbital composition to the 
corresponding bands of Cu2ZnSnS4. 14 In particular, the 
valence band edge has anti-bonding t 2 Cu-d/Se-p* hybrid 
character, with the corresponding bonding states lying in 
the second valence band. Between these regions, at the 
bottom of the first valence band, lies the very dense bun- 
dle of non-bonding e Cu-d orbitals. The two states in the 
lowest conduction band show hybridization among Se-p, 
P-s, and Cu-d orbitals. The charge density distributions 
of these states indicate significant anti-bonding character 
between P and Se ions. This will have important effects 
on relaxation calculations, as discussed below. 

To obtain the optical properties of Cu3PSe4, the 
complex dielectric tensor e has been calculated 15 in 
the random phase approximation from HSE wavefunc- 
tions without including local effects, which are gen- 
erally small. 16 The direction-averaged dielectric func- 
tion Tr(e)/3 and optical absorption coefficient a are 
shown in Figure 3. Optical absorption is compared 
with experimental data 17,18 for GaAs and polycrystallinc 
CuIni_ 2; Ga 2 ;Se2 for x = 0.31. (In photovoltaic applica- 
tions, CIGS stoichiomctries near x = 0.3 yield the most 
efficient power conversion. 19 ) For 520 < A < 660 nm, the 





FIG. 2. Band structure of Cu3PSe4 determined by HSE cal- 
culation. The direct bandgap of 1.38 eV lies at Y. 
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FIG. 1. (Color online) Measured diffuse reflectance of 
Cu3PSe4. Inset: The wurtzite-based enargite structure. For 
Cu 3 PSe 4 , a = 7.685 A, b = 6.656 A, c = 6.377 A. 



FIG. 3. (Color online) Real (solid) and imaginary (long- 
dashed) parts of dielectric function e averaged over principal 
axes. The absorption coefficient a is shown (scale on right) for 
Cu3PSe 4 (red solid), GaAs (brown short-dashed), and CIGS 
(green dashed). 



absorption in Cii3PSe4 is calculated to be higher than in 
GaAs. The greater absorption of CIGS is primarily due 
to its absorption curve being shifted to lower energies by 
its smaller bandgap (1.17 eV). While allowing absorp- 
tion of a greater fraction of the solar spectrum, the CIGS 
bandgap limits its solar cell voltage — and its maximum 
theoretical efficiency — relative to larger bandgap materi- 
als. Absorption is above 5 x 10 4 cm -1 for wavelengths 
A < 630 nm (faj > 2.0 eV) 20 and above 1 x 10 5 cur 1 for 
A < 480 nm (hw > 2.6 eV). 

As discussed previously, the HSE bandgap using the 
experimental atomic and lattice parameters agrees well 
with the measured bandgap. However, using GGA or lo- 
cal density approximation (LDA) relaxed structural pa- 
rameters, we calculate HSE bandgaps approximately 30% 
smaller (TAB. I). While the bandgap is expected to 
depend on the lattice parameters, such a large change 
is unusual, but can readily be understood as follows: 
An HSE ion relaxation leaves the experimental structure 
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TABLE I. Comparison of bandgap and structural quantities 
for different exchange-correlation functionals. Both ion-only 
and full ion+lattice relaxations are performed. All bandgaps 
are predicted to be direct and at T unless noted otherwise. E s 
(HSE) is the HSE-calculated direct bandgap. E g (X) denotes 
the direct bandgap given by measurement or by the functional 
used in the relaxation (LDA, GGA, or HSE). V denotes cell 
volume and the bond lengths are averages. LDA and GGA 
calculations using the unrelaxed experimental structure result 
in bandgaps of 0.21 and 0.29 eV, respectively. 



relaxation 


E g (HSE) 


E s (X) 


V 


(P-Se) 


(Cu-Se) 




(eV) 


(eV) 


(A 3 ) 


(A) 


(A) 


LDA ion 


0.94 


0.08 a 


326.19 


2.30 


2.39 


LDA full 


1.10 b 


0.06 


310.59 


2.27 


2.35 


GGA ion 


1.02 


0.03 


326.19 


2.29 


2.40 


GGA full 


0.90 


0.08 a 


338.25 


2.31 


2.43 


HSE ion 


1.42 


1.42 


326.19 


2.23 


2.42 


exp. 2 


1.38 


1.40 


326.19 


2.24 


2.41 



a An indirect bandgap is predicted (0.04 eV, VBM along T-Y in 
both cases). 

° An indirect bandgap is predicted (1.04 eV, CBM along T-Z). 



and bandgap nearly unchanged. The LDA and GGA 
relaxations, however, result in a lengthening of the P- 
Sc bond, accompanied by large reductions of the HSE 
bandgap. The observed relation between the increase of 
the P-Sc bond length and the underestimate of the cal- 
culated bandgaps is a direct consequence of the signifi- 
cant P-s/Se-p* anti-bonding character of the states at the 
conduction band edge: an underestimate of the bandgap 
(as in LDA/GGA) moves the P-s/Se-p* anti-bonding or- 
bitals to lower energies, resulting in a smaller P-Se bond 
energy and a correspondingly longer bond. A similar 
bond length dependence of the bandgap of CuInSe2 is 
discussed in Rcf. [21]. 

A further structural analysis examines possible al- 
ternatives to the enargite structure used in the pre- 
ceding calculations. While only the enargite form of 
Cu3PSe4 has been observed experimentally, it is possi- 
ble that competing structural phases could form under 
synthesis techniques appropriate for device construction. 
To examine structural stability, we have performed a set 
of GGA calculations (Table II) comparing the heat of 
formation AH for Cu3PSe4 placed in several structures 
which are manifested by other compounds having the 
form A^B 5+ Cl~. The AH values of the famatinite 13 and 
enargite structures differ negligibly. This is not surpris- 
ing, as the essential difference of the structures lies only in 
the arrangement of the (Cu,P)Se4 tetrahedra: zincblende 
for famatinite and wurtzite for enargite. Using the lattice 
and atomic parameters obtained from fully relaxed GGA 
calculations, we have calculated HSE bandgaps of 0.90 
eV (enargite) and 0.84 eV (famatinite) for these poten- 
tially competing structures. Though artificially low due 
to the GGA relaxations, the similarity of these bandgaps 
indicates that the presence of a famatinite phase would 



TABLE II. Calculated zero temperature heat of formation 
(AH) for Cu3PSe4 in several different imaginable structures. 
AH is relative to fully relaxed elements (allotropes: black P 
and 7-Se) and is given per 8-atom formulaic unit. Volume, 
cell shape, and ion positions are relaxed at fixed symmetry. 



structure 



space group 



AH (eV) 



enargite (CU3PS4) 
famatinite (CusSbS4) 
lazarevicite (CU3ASS4) 
Rb 3 PS 4 
Na 3 PS 4 



Pmn2\ (no. 31) 
742m (no. 121) 
P43m (no. 215) 
Pnma (no. 62) 
P42ic (no. 114) 



-1.324 
-1.309 
-1.106 
-0.500 
0.107 



not greatly alter the optical properties of CU3PSC4. 

In conclusion, the measurements and calculations re- 
ported here indicate that Cu3PSe4 has optical properties 
which make it viable for photovoltaic applications. Fur- 
ther measurement and calculation of its properties, ther- 
modynamic stability, and primary defects will be subse- 
quently reported, and will further analyze the potential 
of Cu3PSe4 as a photovoltaic or photoelectronic mate- 
rial. This work has been supported by National Science 
Foundation grant SOLAR DMS-1035513. 
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